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1 Introduction 


The attractor mechanism for supersymmetric (BPS) black holes was discov- 
ered in 1995 [i]: at the horizon of a supersymmetric black hole, the moduli are 
completely determined by the charges of the black hole, independent of their 
asymptotic values. In 2005, Sen showed that all extremal black holes, both su- 
persymmetric and non-supersymmetric (non-BPS), exhibit attractor behavior 
[2]: it is a result of the near-horizon geometry of extremal black holes, rather 
than supersymmetry. Since then, non-BPS attractors have been a very active 
field of research (see for instance [3] Æ [5] [6 [7] 8] O (LO) (11) 03 (13) (74) [T5)). 
In particular, a microstate counting for certain non-BPS black holes was pro- 
posed in [16]. Moreover, a new extension of topological string theory was 
suggested to generalize the Ooguri-Strominger-Vafa (OSV) formula so that it 
also applies to non-supersymmetric black holes [I7]. 

Both BPS and non-BPS attractor points are simply determined as the 
critical points of the black hole potential Vex [I8] [7]. However, it is much 
easier to solve the full BPS attractor flow equations than to solve the non- 
BPS ones: the supersymmetry condition reduces the second-order equations 
of motion to first-order ones. Once the BPS attractor moduli are known in 
terms of D-brane charges, the full BPS attractor flow can be generated via 
a harmonic function procedure, i.e., by replacing the charges in the attractor 
moduli with corresponding harmonic functions: 


tpps(x) =tgpg(p' > H"(x),q1 > Hr(æ)) (1) 


In particular, when the harmonic functions (H? (æ), H;(a)) are multi-centered, 
this procedure generates multi-centered BPS solutions [19]. 

The existence of multi-centered BPS bound states is crucial in under- 
standing the microscopic entropy counting of BPS black holes and the exact 
formulation of OSV formula [20]. One can imagine that a similarly important 
role could be played by multi-centered non-BPS solutions in understanding 
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non-BPS black holes microscopically. However, the multi-centered non-BPS 
attractor solutions have not been constructed until [2I], on which this talk is 
based. In fact, even their existence has been in question. 

In the BPS case, the construction of multi-centered attractor solutions 
is a simple generalization of the full attractor flows of single-centered black 
holes: one needs simply to replace the single-centered harmonic functions in 
a single-centered BPS flow with multi-centered harmonic functions. However, 
the full attractor flow of a generic single-centered non-BPS black hole has not 
been solved analytically, due to the difficulty of solving second-order equations 
of motion. Ceresole et al. obtained an equivalent first-order equation for non- 
BPS attractors in terms of a “fake superpotential,” but the fake superpotential 
can only be explicitly constructed for special charges and asymptotic moduli 
[22] [23]. Similarly, the harmonic function procedure was only shown to apply 
to a special subclass of non-BPS black holes, but has not been proven for 
generic cases [T]. 

In this talk, we will develop a method of constructing generic black hole 
attractor solutions, both BPS and non-BP%S, single-centered as well as multi- 
centered, in a large class of 4D M = 2 supergravities coupled to vector- 
multiplets with cubic prepotentials. The method is applicable to models for 
which the 3D moduli spaces obtained via c*-map are symmetric coset spaces. 
All attractor solutions in such a 3D moduli space can be constructed alge- 
braically in a unified way. Then the 3D attractor solutions are mapped back 
into four dimensions to give 4D extremal black holes. 

The outline of the talk is as follows. Section 2 lays out the framework and 
presents our solution generating procedures; section 3 focuses on the theory 
of 4D N = 2 supergravity coupled to one vector-multiplet, and shows in 
detail how to determine the attractor flow generators; section 4 then uses 
these generators to construct single-centered attractors, both BPS and non- 
BPS, and proves that generic non-BPS solutions cannot be generated via the 
harmonic function procedure; section 5 constructs multi-centered solutions, 
and shows the great contrasts between BPS and non-BPS ones. We end with 
a discussion on various future directions. 


2 Framework 


2.1 3D Moduli Space Msp 


The technique of studying stationary configurations of 4D supergravities by 
dimensionally reducing the 4D theories to 3D non-linear g-models coupled to 
gravity was described in the pioneering work [24]. The 3D moduli space for 
4D N = 2 supergravity coupled to ny vector-multiplets is well-studied, for 
example in [28]. Here we briefly review the essential points. 

The bosonic part of the 4D action is: 
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S=- | dey -g® [R — 2Ggdt Axad — F! AGI (2) 
T 


where I = 0,1...ny, and Gr = (ReN)rzF7 + (ImN) 7 * F7. For a theory 
endowed with a prepotential F(X), Mrz = Fry + 2i E X ee 
Fry = ôrðJ F(X) [28]. We will consider generic stationary solutions, allowing 


non-zero angular momentum. The ansatz for the metric and gauge fields are: 


where 


ds? = —e?” (dt + w)? +67 g dada? (3) 
Al = Ai(dt +w) + A7 (4) 


where ga» is the 3D space metric and bold fonts denote three-dimensional 
fields and operators. The variables are 3ny + 2 scalars {U, t’, t, Aj}, and 
ny +2 vectors {w, A7}. 

The existence of a time-like isometry allows us to reduce the 4D theory to 
a 3D non-linear o-model on this isometry. Dualizing the vectors {w, A7} to 
the scalars {o, By}, and renaming Aj as A‘, we arrive at the 3D Lagrangian, 
which is a non-linear o-model minimally coupled to 3D gravity |] 


1 
C= z V8 7R + 006" Ol" Gmn) (5) 
where @” are the 4(ny + 1) moduli fields {U, t’, t, o, A?, Br}, and gmn is the 
metric of the 3D moduli space M3p, whose line element is: 


ds? = dU? + ye (do + A'dB, — B,dA’)? + gg(t, dt’ dl 
1 -i 
+56 7 [Um (dB; + NixdA®) - (ABs +N jrdA*)] (6) 


The resulting M3p is a para-quaternionic-Kahler manifold, with special 
holonomy Sp(2, R) x Sp(2ny +2, R) [29]. It is the analytical continuation of the 
quaternionic-Kahler manifold with special holonomy U Sp(2, R) x USp(2ny + 
2, R) studied in [26]. Thus the vielbein has two indices (a, A), transforming un- 
der Sp(2,R) and S'p(2ny + 2,R), respectively. The para-quaternionic vielbein 
is the analytical continuation of the quaternionic vielbein computed in [26]. 
This procedure is called the c*-map [29], as it is the analytical continuation 
of the c-map in [25] [26] 

The isometries of the M3p descends from the symmetry of the 4D system. 
In particular, the gauge symmetries in 4D give the shift isometries of M3p, 
whose associated conserved charges are: 


qrdt = Jar = Pg:—B;P,,  p'dr = Jp, = Pp, +A! Po, adr = J, = P, 
(7) 


1 Note that the black hole potential term in 4D breaks down into kinetic terms of 
the 3D moduli, thus there is no potential term for the 3D moduli. 
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where the {P,, Par, Pg,} are the momenta. Here 7 is the affine parameter 
defined as dr = — *3 sin @d0d¢. (p,q ;) are the D-brane charges, and a the 
NUT charge. A non-zero a gives rise to closed time-like curves, so we will set 
a = 0 from now on. 


2.2 Attractor Flow Equations 


The E.O.M. of 3D gravity is Einstein’s equation: 


1 1 
Rab — z8oR = Klan = K(0a™ Ob" Gmn = g 8Bab3c HO H" gmn) (8) 


and the E.O.M. of the 3D moduli are the geodesic equations in M3p: 
VaV G" + Tipa pOH = 0 (9) 


It is not easy to solve a non-linear o-model that couples to gravity. How- 
ever, the theory greatly simplifies when the 3D spatial slice is flat: the dy- 
namics of the moduli are decoupled from that of 3D gravity: 


Tab =0 = 0a oP gmn and 80t” + TR ad OHO (10) 


In particular, a single-centered attractor flow then corresponds to a null 
geodesics in Mp: ds? = døde” gmn = 0. 

The condition of the 3D spatial slice being flat is guaranteed for BPS 
attractors, both single-centered and multi-centered, by supersymmetry. Fur- 
thermore, for single-centered attractors, both BPS and non-BPS, extremality 
condition ensures the flatness of the 3D spatial slice. In this paper, we will 
impose this flat 3D spatial slice condition on all multi-centered non-BPS at- 
tractors we are looking for. They correspond to the multi-centered solutions 
that are directly “assembled” by single-centered attractors, and have prop- 
erties similar to their single-centered constituents: they live in certain null 
totally geodesic sub-manifolds of M3p. We will discuss the relaxation of this 
condition at the end of the paper. 

To summarize, the problem of finding 4D single-centered black hole attrac- 
tors can be translated into finding appropriate null geodesics in M3p, and 
that of finding 4D multi-centered black hole bound states into finding cor- 
responding 3D multi-centered solutions living in certain null totally geodesic 
sub-manifold of M3p. 

The null geodesic that corresponds to a 4D black hole attractor is one that 
terminates at a point on the U — —oo boundary and in the interior region 
with respect to all other coordinates of the moduli space M3p. However, it 
is difficult to find such geodesics since a generic null geodesic flows to the 
boundary of Msp. For BPS attractors, the termination of the null geodesic 
at its attractor point is guaranteed by the constraints imposed by supersym- 
metry. For non-BPS attractor, one need to find the constraints without the 
aid of supersymmetry. We will show that this can be done for models with 
Msp that are symmetric coset spaces. Moreover, the method can be easily 
generalized to find the multi-centered non-BPS attractor solutions. 
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2.3 Models with Msp Being Symmetric Coset Spaces 


A homogeneous space M is a manifold on which its isometry group G acts 
transitively. It is isomorphic to the coset space G/H, with G being the isom- 
etry group and H the isotropy group. For M3p = G/H, H is the maximal 
compact subgroup of G when one compactifies on a spatial isometry down to 
(1,2) space, or the analytical continuation of the maximal compact subgroup 
when one compactifies on the time isometry down to (0,3) space. 

The Lie algebra g has Cartan decomposition: g = h @k where 


(h,h] =h [h,k] =k (11) 
When G is semi-simple, the coset space G/H is symmetric, meaning: 
({k,k] =h (12) 


The building block of the non-linear o-model with symmetric coset space 
Msp as target space is the coset representative M, from which the left- 
invariant current is constructed: 


J=M—dM=A+dn (13) 


where J, is the projection of J onto the coset algebra k. The lagrangian 
density of the o-model with target space G/H is then given by Jy as: 


L= Tr( J A *3 Jk) (14) 


The symmetric coset space has the nice property that its geodesics M (7) 
are simply generated by exponentiation of the coset algebra k: 


M(t) = Moe"? with kek (15) 


where Mo parameterizes the initial point of the geodesic, and the factor 4 in 
the exponent is for later convenience. A null geodesic corresponds to |k|? = 0. 
Therefore, in the symmetric coset space M3p, the problem of finding the 
null geodesics that terminate at attractor points is translated into finding the 
appropriate constraints on the null elements of the coset algebra k. 

The theories with 3D moduli spaces M3p that are symmetric coset spaces 
include: D-dimensional gravity toroidally compactified to four dimensions, all 
4D N > 2 extended supergravities, and certain 4D M = 2 supergravities 
coupled to vector-multiplets with cubic prepotentials. The entropies in the 
last two classes are U-duality invariant. In this talk, we will focus on the last 
class. The discussion on the first class can be found in [21]. 


Parametrization of Msp 


The symmetric coset space M3p = G/H can be parameterized by exponenti- 
ation of the solvable subalgebra solv of g: 
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The solvable subalgebra solv is determined via Iwasawa decomposition of g. 
Being semi-simple, g has Iwasawa decomposition: g = h@a@n, where a 
is the maximal abelian subspace of k, and n the nilpotent subspace of the 
positive root space X+ of a. The solvable subalgebra solv = a@n. Each point 
ġo” in M3p corresponds to a solvable element ©(¢) = e°% , thus the solvable 
elements can serve as coset representatives. 

We briefly explain how to extract the values of moduli from the coset 
representative M. Since M is defined up to the action of the isotropy group 
H, we need to construct from M an entity that encodes the values of moduli 
in an H-independent way. The symmetric matrix S defined as: 


S = MS MT (17) 


has such a property, where Sọ is the signature matrix] Moreover, as the 
isometry group G acts transitively on the space of matrices with signature So, 
the space of possible S is the same as the symmetric coset space M3p = G/H. 
Therefore, we can read off the values of moduli from S' in an H-independent 
way. 

The non-linear o-model with target space M3p can also be described in 
terms of S instead of M. First, the left-invariant current of S is Jg = ST tds, 
which is related to Jẹ by: 


Js = S71dS = 2(So MT)! Jkl So MT) (18) 


The lagrangian density in terms of S' is thus £ = {Tr(Js^*3Jg). The equation 
of motion is the conservation of current: 


V-J=V-(S~'VS) =0 (19) 


where we have dropped the subscript S in Js, since we will only be dealing 
with this current from now on. 


2.4 Example: ny = 1 


In this talk, we will perform the explicit computation only for the simplest 
case: 4D N = 2 supergravity coupled to one vector-multiplet. The generaliza- 
tion to generic ny is straightforward. The 3D moduli space M3p for ny = 1 
is an eight-dimensional quaternionic kahler manifold, with special holonomy 
Sp(2, R) x Sp(4, R). Computing the killing symmetries of the metric (6p with 
ny = 1 shows that it is a coset space Go(2)/(SL(2, R) x SL(2,R))P} Figure [1] 
shows the root diagram of Gai) in its Cartan decomposition. The six roots on 


2 In all systems considered in the present work, the isotropy group H is the maximal 
orthogonal subgroup of G: HSo HT = So, for any H € H. Therefore, S is invariant 
under the H-action M — MH. 

3 Other work on this coset space has appeared recently, including [33]. 
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Fig. 1. Root Diagram of G'g/2) in Cartan Decomposition. 


the horizontal and vertical axes {L7 , L3, LF, L3} generate the isotropy sub- 
group H = SL(2, R)» x SL(2, R)». The two vertical columns of eight roots aa 
generate the coset algebra k, with index a labeling a spin-1/2 representation 
of SL(2, R), and index A a spin-3/2 representation of SL(2, R)». 

The Iwasawa decomposition, g = h @ solv with solu = a È n, is 
shown in Figure |2| The two Cartan generators {u,y} form a, while n is 
spanned by {x,o0,A°,A',Bi,Bo}. {u, y} generates the rescaling of {u, y}, 
where u = e74, and {x,o0,A°,A',B,,Bo} generates the translation of 
{x, o, A°, A’, B,, Bo} [27]. 

The moduli space M3p can be parameterized by solvable elements: 


Xlo) = elinu)u/2+(Iny)y pext+AT A+B Bitoo (20) 


The symmetric matrix S can then be expressed in terms of the eight moduli 
on: 

S($) = EA) E)" (21) 
which shows how to extract the values of moduli from S' even when S is not 
constructed from the solvable elements, since it is invariant under H-action. 
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Fig. 2. Root Diagram of Gə) in  Iwasawa Decomposition. 


{u, y, x, øo, A?, A!, B1, Bo} generates the solvable subgroup Solv. 


3 Generators of Attractor Flows 


In this section, we will solve 3D attractor flow generators k as in (15). We 
will prove that extremality condition ensures that they are nilpotent elements 
of the coset algebra k. In particular, for ny = 1, both BPS and non-BPS 
generators are third-degree nilpotent. However, despite this common feature, 
kgps and kyp differ in many aspects. 


3.1 Construction of Attractor Flow Generators 
Construction of kgps 


Since the 4D BPS attractor solutions are already known, one can easily obtain 
the BPS flow generator kgpg in the 3D moduli space M3p. 

The generator kgps can be expanded by coset elements aaa as kppg = 
Qo AOMA, where C°“ are conserved along the flow. On the other hand, since 
the conserved currents in the homogeneous space are constructed by project- 
ing the one-form valued Lie algebra g7™t-dg onto k, a procedure that also gives 
the vielbein: Jk = g~'dg|k = adaaV, the vielbein V®^4 are also conserved 


along the flow: 4 (veer) = 0. Since both the expansion coefficients C°4 
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and the vielbein V°4 transform as (2,4) of SL(2, R)» x (SL(2,R), and are 
conserved along the flow, they are related by: 


C (22) 


up to an overall scaling factor. 

In terms of the vielbein V°“, the supersymmetry condition that gives the 
BPS attractors is: V°4 = z°V4 [29] B0 B3]. Using (22), we conclude that the 
3D BPS flow generator kgps has the expansion 


kpps = dqa2z*C4 (23) 


A 4D supersymmetric black hole is labeled by four D-brane charges 
(p°, pt, qi, qo). A 3D attractor flow generator kgps has five parameters {O4, z}. 
As will be shown later, z drops off in the final solutions of BPS attractor flows, 
under the zero NUT charge condition. Thus the geodesics generated by kpps 
are indeed in a four-parameter family. 

kgps can be obtained by a twisting procedure as follows. First, define 
a k9 ps Which is spanned by the four coset generators with positive charges 
under S'L(2,R)p: 

k% pg = a404 (24) 


then, conjugate k9 pg With lowering operator L;, : 
kpps =e "nk ppgevn (25) 
Using properties of ke pg, it is easy to check that kgpg is null: 
|keps|? =0 (26) 
More importantly, kgpg is found to be third-degree nilpotent: 
kbps =0 (27) 


A natural question then arises: Is the nilpotency condition of kgps a result of 
supersymmetry or extremality? If latter, we can use the nilpotency condition 
as a constraint to solve for the non-BPS attractor generators kyg. We will 
prove that this is indeed the case. 


Extremality implies nilpotency of flow generators 


We will now prove that all attractor flow generators, both BPS and Non-BPS, 
are nilpotent elements in the coset algebra k. It is a result of the near-horizon 
geometry of extremal black holes. 

The near-horizon geometry of a 4D attractor is AdS2 x S?, i.e. 


ee! — /Veu\|sT as TOO (28) 
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As the flow goes to the near-horizon, i.e., as u = e?” — 0, the solvable element 


M = eln u)u/2+ is a polynomial function of 7: 
M(t) ~ ut? ~ rf (29) 


where —£ is the lowest eigenvalue of u. 

On the other hand, since the geodesic flow is generated by k € k via 
M(r) = Moe*’/?, M(r) is an exponential function of 7. To reconcile the two 
statements, the attractor flow generator k must be nilpotent: 


kt = 9 (30) 


where the value of depends on the particular moduli space under consid- 
eration. In Ga2)/SL(2,R)?, by looking at the weights of the fundamental 
representation, we see that Z = 2, thus 


k3 =0 (31) 


The nilpotency condition of the flow generators also automatically guar- 
antees that they are null: 


k’ =0 = rap = Tr(k?) =0 (32) 


Construction of kyng 


To construct non-BPS attractor flows, one needs to find third-degree nilpotent 
elements in the coset algebra k that are distinct from the BPS ones. In the 
real Ga) /SL(2,R)?, there are two third-degree nilpotent orbits in total [85]. 


We have shown that kgpg = e~*4nk9,pce*%», with k} pg spanned by the 
four generators with positive charge under SL(2, R)». 

Since there are only two SL(2,R)’s inside H, a natural guess for ky p is 
that it can be constructed by the same twisting procedure with SL(2, R)» 
replaced by SL(2, R)»: 


kyng =e kl peo with «= hp = Aaa", a,a=1,2 (33) 


where k9p is spanned by the four generators with positive charge under 
SL(2, R)». 

Using properties of k9, p, one can easily show that kyng defined above is 
indeed third-degree nilpotent: 


kłp=0 (34) 


That is, kyp defined in generates non-BPS attractor flows in M3p. 

A 4D non-BPS extremal black hole is labeled by four D-brane charges 
(p°, p', q1, qo). Similar to the BPS case, the 3D attractor flow generator ky g 
has five parameters {C'°’,z}. As will be shown later, z can be determined 
in terms of {C°*} using the zero NUT charge condition, thus the geodesics 
generated by kyp are also in a four-parameter family. 
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3.2 Properties of Attractor Flow Generators 


We choose the representation of Gg(2) group to be the symmetric 7 x 7 
matrices that preserve a non-degenerate three-form Wijk such that nis = 
WijkWstuWmnoe evr” is a metric with signature (4,3) and normalized as 
n? = 1. We decompose 7 as 36341 of SL(3,R) and choose the non-zero 
components of w,3A3A3,3A3A3 and 3893 1, as 


1 
v2 


which gives 7 = dxzady — dz?. Written explicitly, an element of Go(2) Lie 
algebra is 


w = dx, A dzz ^ dx3 + dy! A dy? A dy? — —dzra A dy? A dz (35) 


AÑ  enjakt V2wi, 
g= eth ayy =A; —/2v2 (36) 
— 2v: V2w;, 0 


Here A is a traceless 3x3 matrix. The signature matrix Sp is thus Diag[1, —1,—1,1,—1,—1, 1]. 
The real G2(2) group has two third-degree nilpotent orbits. In both orbits, 
k? is of rank two and has Jordan form with two blocks of size 3. Thus k? can 
be written as 
k? = 5 VaV Cab S0 (37) 


a,b=1,2 
with va null and orthogonal to each other: vg - vp = vT Soup = 0, and Cab 
depends on the particular choice of k. Therefore, k can be expressed as: 


k= 5 (vawe + wave) So (38) 


a=1,2 


where each wa is orthogonal to both vg: Wa ` vp = 0, and Wa satisfies we - 
Wp = Cap. Next we solve for vg and Wa for kgpg and kyg and compare their 
properties. 


Properties of kgps 


The null space of k? is five-dimensional, with va spanning its two-dimensional 
complement. For kgps, goes and were in are solved in terms of C^ 
and z. 

In basis (36), from inspection of k% p s, we find that vBPS can always be 


chosen to have the form] 


vPPS = (Vi, =m V1,0) vB PS = (—V2, MV2, V2) (39) 


“There are some freedom on the choice of (va, Wa): a rotational freedom: 
(Ua, Wa) — (Rabvo, Raw) with R orthogonal; and a rescaling freedom: 
(Va, Wa) — (Var, Wa fr). 
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where nı is a 3D signature matrix 7 = Diag|1, —1,—1], and V, are two three- 
vectors satisfying 
Y-V, =0 Vv- V2 =0 V- V = —1 (40) 
We drop the superscript “BPS” for V, here since, as will be shown later, v? 
can also be written in terms of V,, though in a slightly different form. Note 
that for kgpg, V2 is defined up to a shift of Vj: V2 > Və — cV}, since any linear 
combination of v2? forms a new set of vB?%, 
Written in twistor representation [?] Va are given by the twistors z and u as 


VEP = 22%2P Ve? = z%b + zbut (41) 


where we have used the rescaling freedom to set ztu? — z2u! = 1. Note that 
for kgps, the twistor u is arbitrary, due to the shift freedom of V3. 


The condition w??S . vB PS = 0 dictates that w??% has the form: 


wP? = (WPS, mWPPF,0) wgs = (WPS, mWyP*,0) (42) 
with WPPS solved as: 
(WPPS, ee tad — (PPVU, PoP x.) (43) 
where the totally symmetric P®° is defined in terms of C4 as 
pil = Cc! pile = C? pi22 = C3 p??? = C4 (44) 


BPS span a one-dimensional space (since u is arbitrary) and 


span a four-dimensional space. 


In summary, v 
wBPS 


3.3 Properties of kNB 


(vP wN®) are solved in terms of {C°*, z}. The forms of v are only slightly 


different from those of vP?S: 


uB = (Vi, m Vi, 0) ug” = (V2, =m Va, V2) (45) 


where V, are the same three-vectors given in (41), with one major difference: 
the twistor u is no longer arbitrary, but is determined by C'°“ as 
2 22 
u C 
u C 
5 With the inner product of two three-vectors defined as Va - Ve = VZ m Ve, the 
twistor representation of a three-vector V = (x,y,z) can be chosen as 


Ov = £00 + yo3 + ZO, = TEY Z 
z @g-y 
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since the V2 in v no longer has the shift freedom. 


The forms of wN? are also only slightly different from the BPS ones (42): 
wy? a (WDP, —mWy*, 0) wy? = (WH, mW; ®,0) (47) 
with WN solved in terms of {C°, z, u} as: 


(wre = uĉuf + (Cy? — C1u!)z”zf (48) 
(We? 8 = ea? + u®2f) + (0?! Ce Bi ete (49) 


Since the value of u imposes an extra constraint on the vectors wN B via (46), 
wN® span a three-dimensional space instead of a four-dimensional one as in 
the BPS case (43). In summary, in contrast to the BPS case, v? span a 


a 
two-dimensional space and wÑ? span a three-dimensional one. 


A Single-centered Attractor Flows 


Having solved the attractor flow generators for both BPS and non-BPS case, 
we are ready to construct single-centered attractor flows. A geodesic starting 
from arbitrary asymptotic moduli is given by M(r) = Moe*"/?, which gives 
the flow of S as S(T) = Moe? SoM@,, which in turn can be written as S(T) = 
e*() So, where K(r) is a matrix function. From now on, we use capital K 
to denote the matrix function which we exponentiate to generate attractor 
solutions. 
The current of S is 


1 
J= SVS = S (v+ [VK, K] + 5[VK,K],K] ++- +) So (50) 
The equation of motion is the conservation of currents: V - (S~!VS) = 0, 
which is solved by K(r) being harmonic: 
VE) =0 = K(r)=kr+g (51) 


g parameterizes the asymptotic moduli. Using the H-action, we can adjust 
g such that g € k, and g has the same properties as the flow generator k, 
namely, g? = 0 and g? is of rank two. Therefore, for single-centered flow 
given by S(T) = e*() So, the harmonic matrix function K(r) has the same 
properties as the flow generator k: 


K3(r) =0 and K?(r) rank two (52) 


To find the harmonic K(r) that satisfies the constraints (52), recall that 
the constraints dictate K (rT) to have the form: 


K(r) = 5 (va(T)wa(T)? + wa(T)V¥a(T)") So (53) 
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with vq(7) being null and wa(T) orthogonal to va(T) for all r. Then the con- 
straints (52) can simply be solved by choosing va(T) to be the constant null 
vectors ValT) = Va and wa(T) to be harmonic vectors which are everywhere 
orthogonal to vp: 


LS 
— 


Wa(T) = WaT + Ma with Wa * Ub = Ma * Vb = 0 (54) 


The two 7-vectors w,’s contain the information of the black hole charges, and 
the two 7-vectors m,’s contain that of asymptotic moduli. 

To summarize, the single-centered attractor flow starting from an arbitrary 
asymptotic moduli is generated by S(T) = e* 7) So, with harmonic matrix 
function K(T) = kr + g where 


k= 5 [vaw] + wave |So and g= 5 [vam] + mav?]So (55) 


a=1,2 a=1,2 


Since k and g share the same set of null vectors vg and both wa and mg, are 
orthogonal to vp, g has the same form as that of flow generator k, namely: 


JBPS = Gea2z"G* gnB = € 7" (aaa) e7 (56) 


which guarantees that g is also third-degree nilpotent. Moreover, that g and 
k have the same form implies [[k, g], g] = 0, thus the current is reduced to 


_ So(k + alk, 9) So .. 


J 2 


= (57) 
from which we can solve vg and wa in terms of charges and asymptotic moduli. 

Now that we are able to construct arbitrary attractor flows in the 3D 
moduli space, we can lift them to the 4D black hole attractor solution. First, 
in representation given by (36), the 4D moduli t = x + iy can be extracted 
from the symmetric matrix S via: 


535(7) ole — S35(T)? 


noss S33(T) = S33 (7)? 


(58) 


1 
u = 
V S33(T)S55(T) — S35(7)? 
Since both k and g are third-degree nilpotent, S(T) is a quadratic function of 
T. Moreover, since g has the same form as k, S(T) is composed of harmonic 
functions of r: H4(r) = C47 +Gé4 for BPS attractors and H°*(r) = C%r+ 
G°* for non-BPS attractors [E] Generic single-centered attractor flows with 
arbitrary charges and asymptotic moduli can thus be generated. The attractor 


(59) 


6 Space prohibits listing the rather lengthy result of S(T), readers can consult eq. 
(6.2) and eq. (6.3) of for BPS attractors, and eq. (6.28) for non-BPS ones. 
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moduli are read off from S(T) with T — oo, and asymptotic moduli with 
T > 0. 

The D-brane charges can be read off from the charge matrix defined as 
Q= x fN - J. The 4D gauge currents sit in the current J = S~'VS as: 


a 2 2 
(J31, Js1, J72, J12, J32) = (V2JB,, —V2JB,, gar V2Jao, —2Jo) (60) 


Therefore Q relates to the D-brane charge (p°, p!,q1,qo) and the vanishing 
NUT charge a by 


(Qs1, Qs1, Qr2, Qi) = (VPP, -V2p", 51, V240) Qs = -2a =0 (61) 


4.1 Single-Centered BPS Attractor Flows 


As an example, a single-centered BPS black hole constructed by lifting the 
attractor solution in Mgp is shown in Figure It has D-brane charges 
(pe, pt, qı, qo) = (5,2,7, —3). The four flows, starting from different asymp- 
totic moduli, terminate at the attractor point (tpg, ypzpg) with different 
tangent directions. The reason is that the mass matrix of the black hole po- 
tential Vey at the BPS critical point has two identical eigenvalues, thus there 
is no preferred direction for the geodesics to flow to the attractor point. 

We now discuss in detail how to determine kgpg and ggpg for given 
charges and asymptotic moduli. There are nine parameters in kgps and ggps: 
{C4, G4, z}, since the twistor u is arbitrary. On the other hand, there are 
eight constraints in a given attractor flow: four D-brane charges (p/,q7), the 
vanishing NUT charge a, and the asymptotic moduli (x0, Yo, to) {"] We will 
use these eight constraints to fix C4 and G4 in kgpg and ggps, leaving the 
twistor z unfixed. 

Integrating the current for BPS case produces five coupled equations: 


1 
Qeps = So(keps + z lkeps, 9BPs])So (62) 


where [kpps,gpps| = (C,G)O, with (C, G) = C1G*—3C7G?+3C°G?-C4G", 
and O = — rae en Le. 

In order to show that the BPS flow can be expressed in terms of harmonic 
functions: H(r) = Qr +h, with Q = (p!,qr) and h = (h', hr), we will solve 
geps in terms of h instead of (xo, yo, uo). h relates to the asymptotic moduli 
by 


(£0, Yo, Uo) BPS = (x, Y, u)bps(Q = h) (63) 


and there is one extra degree of freedom to be fixed later. 


T The asymptotic value of u can be fixed to an arbitrary value by a rescaling of 
time and the radial distance. We will set uo = 1. 
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X4,Y4) 


L i l i 1 Ì ıı ı l x 


-1.0 -0.5 0.5 1.0 1.5 2.0 


Fig. 3. BPS flow with charge (p°,p',qi,qo) = (5,2,7,—3) and attractor point 
(a*,y*) = (0.329787, 0.788503). The initial points of each flow are given by (xı = 
1.5, yı = 0.5), (xo 2, Y2 4), (x3 0.2, y3 = 0.1), (x4 = —1, y4 = 3). 


First, for later convenience, we separate from ggps a piece that has the 
same dependence on (h, z) as kgpg on (Q, z): 


IBPS = JBPS,h + A with JBPS,h = kgps(Q > h, z) (64) 


that is, gBPS,h = doaz°G# with GA = C4(Q — h). We can use the unfixed 
degree of freedom in h to set (C, Gn) = 0, so that (62) simplifies into 


1 
Qzaps = So(kaps + glkeps: A]) So (65) 
A can then be determined using the three constraints from (63) and the zero 
NUT charge condition in (65): A = aaz“ E^ with E! = -F3 = =z and 
E? = -Et = Gz The form of A will ensure that the twistor z drops off in 


the final attractor flow solution written in terms of Q and A. 
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The remaining four conditions in the coupled equations determine C4 
as functions of D-brane charges and the twistor z: C4 = C4(Q, 2) P|T hen GA 
are given by G = C4(Q — h, z). The product (C4, G#) is proportional to 
the symplectic product of (p’,qr) and (h/, hr): 


(04, GA) = ae <Q,h> where < Q, h >= p°hotpthi—qih'—qoh® 
(66) 
The condition (C4,GA) = 0 is then the integrability condition on h: < 
Q,h >=0. 
BPS attractor flows in terms of (p7, qr) and (h?, hr) are obtained by substi- 
tuting solutions of C4 (Q, z) and GÅ (h, z) into the flow of S(T). The attractor 
moduli are determined by the charges as: 


0 1qı O ml AL 
+ Ey J 3 2°30 
Paty 4 J Ja(p", pt, &, go) (67) 


CAPR "DRY + pee] NRPS Tp) + poe] 


where J4(p°, p', 91, qo) is the quartic E7(7) invariant: 


Ja(p°, p*, qi, qo) = 3(p'q1)?—-6(p° Go) (p' a1) —(P° G0)? —4(p") go +4" (q1)? (68) 


thus J4(p°,p',4,qo) is the discriminant of charges. Charges with positive 
(negative) J4(p°,p',4,qo) form a BPS (non-BPS) black hole. The attractor 


value of u is ugpg = 1/,/Ja(p°, pt, 4, qo). The constraint on h from uo = 1 is 
then J4(h°, ht, 4+, ho) = 1. The attractor moduli match those from Type 
II string compactified on diagonal Tî, with q, > F. 

Now we will prove that the BPS attractor flows constructed above can 
indeed be generated by the “naive” harmonic function procedure, namely, by 
replacing charges Q in the attractor moduli with the corresponding harmonic 
functions Qr + h. First, using the properties of A, the flow of t = x + iy can 
be generated from the attractor moduli by replacing kgpg with the harmonic 


function kgpsT + 9BPS,nh: 


= teps(T) = thps(kBPs — kr + 9BPS,h) (69) 


Then, since kgps and ggps,n share the same twistor z, this is equivalent to 
replacing C4 with harmonic functions H4(7) = C47 + GA while keeping the 
twistor z fixed: 

ipps(T) =tepg(C4 > OTF Gs z) (70) 


Finally, since C^ is linear in Q and GA linear in h, and since z drops off after 
plugging in the solutions C4(Q, z) and G#(h, z), we conclude that the flow of 
tgps(T) is given by replacing the charges Q in the attractor moduli with the 
corresponding harmonic functions Qr + h: 


tpps(T) = teps(Q > QT +h) (71) 
8 See eq. (6.18) of for the full solutions. 
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4.2 Single-Centered Non-BPS Attractor Flows 


A non-BPS attractor flow with generic charges and asymptotic moduli can 
be generated using the method detailed earlier. Figure [4] shows an example 
of non-BPS attractor flow with charges (p°, p', q1, qo) = (5, 2,7,3). Note that 
J4(5,2,7/3, 3) < 0, so this is indeed a non-BPS black hole. 


(x1,y1) 


(x4, Y4) 
L (x2, y2) 
(x5,y5) L 
(X3 ya) i] — l x 
—0.5 0.5 1.0 1.5 


Fig. 4. Non-BPS flow with charges (p°, p', qi, qo) = (5,2,7,3) and attractor point 
(x*,y*) = (—0.323385, 0.580375). The initial points of each flow are given by: (xı = 
0.539624, yı = 5.461135), (z2 = 1.67984, y2 = 0.518725), (zs = —0.432811, ys = 
0.289493), (va = 1.28447, ya = 1.49815), (ws = —0.499491, ys = 0.181744). 


Unlike the BPS attractor flows, all non-BPS flows starting from different 
asymptotic moduli reach the attractor point with the same tangent direction. 
The reason is, unlike the BPS case, the mass matrix of the black-hole potential 
Vex at a non-BPS critical point has two different eigenvalues. The common 
tangent direction for the non-BPS flows corresponds to the eigenvector asso- 
ciated with the smaller mass. 
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Now we discuss how to determine ky g and gyp for given D-brane charges 
and asymptotic moduli. Unlike the BPS case, there are only eight parameters 
in ky p and gyp: the two twistors {z, u} and {C°*, G°*} under the constraints 
u= oit: = Gi. On the other hand, there are still eight constraints in a given 
non-BPS attractor flow as in the BPS case. Therefore, while kgps and ggps 
can parameterize black holes with arbitrary (p!, qr) and (£o, yo) while leaving 
{z,u} free, all the parameters in ky p and gy zp, including {z, u}, will be fixed. 

Another major difference from the BPS case is that 


[kn B,gNnB] = 0 (72) 


guaranteed by the form of vP and wN®. Thus the charge equation 
becomes simply 
Que = So(knz) So (73) 


Unlike the BPS case, gy g does not enter the charge equations, thus cannot 
be used to eliminate the dependence on the twistor z. The three degrees 
of freedom in gyp are simply fixed by the asymptotic moduli (zo, yo) and 
uo = 1, without invoking the zero NUT charge condition. The four D-brane 
charges equations in determine C°* = c7(Q, 2) FI which then fixes u 
via u = Siz. Finally, the zero NUT charge condition imposes a degree-six 


equation on twistor z: 
p°2°+6p' 2° — (3p? +4q1)z*—4(3p*—2q0)2?+ (3p°+4q1)z7+6p'z—p° =0 (74) 


Similar to the BPS case, the full non-BPS attractor flow can be gener- 
ated from the attractor moduli by replacing C°* with the harmonic function 
H°” (T) = C%'r + G°*, while keeping z fixed as in (70): 


tNB(T) = typ(C™ — Cr + G°", z) (75) 


However, there are two important differences. First, the harmonic functions 
H°* have to satisfy the constraint {™] 


H??? (T) C22 G22 
= u = 


H2(r) Cl ` G2 


(76) 


Second, unlike the BPS flow, a generic non-BPS flow cannot be given by 
the “naive” harmonic function procedure: 


tna(tT) Æ tNB(Q > QT +h) (77) 


? See eq. (6.35) of for the full solution. 

10 This does not impose any constraint on the allowed asymptotic moduli since there 
are still three degrees of freedom in G** to account for (xo, yo, uo). We will see 
later that its multi-centered counterpart helps impose a stringent constraint on 
the allowed D-brane charges in multi-centered non-BPS solutions. 
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The reason is that the twistor z in a non-BPS solution is no longer free as in 
the BPS case, but is determined in terms of D-brane charges via (74). Thus 
replacing Q with QT+h, for generic Q and h, would not leave z invariant. That 
is, replacing C° in the attractor moduli with harmonic functions H®°(r) is 
not equivalent to replacing the charges Q with H = Qr + h as in the BPS 
case (71). 

It is interesting to find the subset of non-BPS single-centered flows that 
can be constructed via the “naive” harmonic function procedure. The ny = 1 
system can be considered as the STU model with the three moduli (S, T, U) 
identified. Since the STU model has an SL(2,Z)? duality symmetry at the 
level of E.O.M., the ny = 1 system has an SL(2, Z) duality symmetry coming 


i ia 7 a (ab ab ab 
from identifying these three SL(2,Z)’s, namely, I’ = p AGE AGE a) 


with ad — bc = 1. The modulus t = x + iy transforms as t — I't = 
the transformation on the charges is given by [34]. 

Given an arbitrary charge Q, there exists a transformation Î Q such that 
Q= ÎoQao for some D4-D0 charge system Q40 = (0, pt, 0, qo). The solution 
of with charge Q = Ê QQao has a root z = ee vanie independent of Q40. 
Thus for arbitrary h4o = (0, ht,0, ho), replacing Q with Qr + Îohaso would 
leave the twistor z invariant. We thus conclude that the non-BPS single- 
centered attractor flows that can be generated from their attractor moduli 
via the “naive” harmonic function procedure are only those with (Q, h) being 
the image of a single transformation Î on a D4-D0 system (Quo, hao): 


at+b 
ct+d’ 


and 


twa(t) = typ(L'Qa0 > PQaot + Pho) (78) 


5 Multi-Centered Attractor Flows 


Similar to the single-centered attractor solutions, the multi-centered ones are 
constructed by exponentiating harmonic matrix functions K (a): 


S(x) = e*) So (79) 


Recall that for single-centered attractors, using the H-action on g, K(r) = 
kr + g can be adjusted to have the same properties as the flow generator k 
as in (52). For BPS multi-centered solutions, supersymmetry guarantees that 
the matrix function K(x) also has the same properties as the generator k: 


K(x) =0 and K?(a) rank two (80) 


We will impose these constraints on all non-BPS multi-centered solutions as 
well, since presently we are more interested in the multi-centered solutions that 
are “assembled” by individual single-centered attractors and thus have similar 
properties to their single-centered constituents. It is certainly interesting to 
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see if there exist non-BPS multi-centered solutions with K(x) not sharing the 
constraints satisfied by the flow generator ky gp. 

The harmonic matrix function K(a) satisfying all the above constraints is 
solved to be: 


K(x) = 3 = +g (81) 


ky = 5 [valWa) E +(wa)iv?]So and g= 5 [vam] +mav?]So (82) 
a=1,2 


a=1,2 


with vg being the same two constant null vectors in single-centered k, and the 
7-vectors (Wa); contain the information of the D-brane charges of center-i, 
and the two 7-vectors m,’s contain that of asymptotic moduli. Both (wa)i 
and Ma are orthogonal to vp. Since vg only depends on the twistor {z,u}, 
and wa are linear in C4 or C, the above generating procedure is equivalent 
to replacing C4 and C°®* with the multi-centered harmonic functions H4(a) 
and H°“(a) while keeping the twistor {z,u} fixed. 


5.1 Multi-Centered BPS Attractors 


Using Q; to denote the charge matrix of center-i, we have 5N coupled equa- 
tions from Q; = x JV J: 


KBPS,i, KBPS,j 
[k BPS, BPS, g 


1 1 
nF | keps t kupia 83 
QBpPs, olkBps,i + slkepsi. gers) + 5 2 Beca] o (83) 


We now show in detail how to determine k gps, and ggps for given charges 
and asymptotic moduli using equation (83). There are 4(N+1)+1 parameters 
in kpps,; and ggps: Gs GA Zy, since the twistor u is arbitrary. Different 
from the single-centered BPS case, there are also 3N — 3 degrees of freedom 
from the positions of centers on L.H.S. of (83). On the other hand, there 
are 5N +3 constraints in a given BPS multi-centered attractor: 4N D-brane 
charges (p!,qr,;), N vanishing NUT charges, the asymptotic moduli (xo, yo) 
and uo = 1. We will use these 5N +3 constraints to fix the 4(N +1) parameters 
{CA, G4} in kgps,i and ggps, and impose N — 1 constraints on the distances 
between the N centers, while leaving the twistor z free. 

First, integrating V - J over the sphere at the infinity gives the sum of 
the above N matrix equations: Qb s = So(kfbs + $[k8bs,9BPs]) So, which 
is the same as the charge equation for a single-center attractor with charge 

tt g. This determines g to be g = gn + A, same as the single-centered case 
as in (64), using the three asymptotic moduli (£o, yo, uo) and the constraint 
of zero total NUT charge. The h’s are fixed by the asymptotic moduli and the 
integrability condition < Qh g, h >= 0. 
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It is easy to see that the solutions of CA are simply given by the single- 
centered solutions C4 = C4(Q,z) with Q replaced by Q;. Thus the flow 
generator of each center kgps, (given by kgps; = Gear) satisfies 


1 
Qeps,i = Sol(keps,i + al BPs,i, Al) So (84) 


which is the multi-centered generalization of the single-centered condition 


(65). 


Using the solutions of kgps,, and ggps, the charge equations (83) become 


< Qepsi,QpBps,j > 
|x; — æj| 


(85) 


)O | So 


1 
Qeps,i = So | KBPS, + glkeps.i: A) +(< Qzpsi,h>+> 
j 
from which we subtract (84) to produce the integrability condition 


< QBPS,i, QBPS,j > 
|x; — j| 


< Qppsi,h > +5 =0 (86) 
j 


a l The sum of the N equations in the integrability condition reproduces 
= the constraint on h: < Qb, h >= 0. Thus the remaining N — 1 equations 
impose N — 1 constraints on the relative positions between the N centers. The 


. i k . 
angular momentum J, defined via w; = 2€;;,J7 fs as r — oo, is non-zero: 


1 Ti — Lj 
J= 5 > Jay wa] BPS QBPs,;) (87) 


AI Thus we have shown that our multi-centered BPS attractor solutions repro- 
nu duce those found in [I9]. Same arguments as in the single-centered BPS case 
shows that multi-centered BPS attractors can be generated by replacing the 
= charges in the attractor moduli with corresponding multi-centered harmonic 
T functions: 
= tpps() = thps(QBPS > 5 gers: + h) (88) 
a 


z- gzil 


5.2 Multi-Centered Non-BPS Attractors 


A multi-centered non-BPS attractor has 3(N+1)+2 parameters inside its non- 
BPS generators {knB i, gNB}: {Ce%,G°*} under the constraint plus two 
twistors {z,u}. Given {kyNB i, gyng} in terms of {09°, G°*, z, u}, the non-BPS 
multi-centered solution is the same as the single-centered one with H°*(r) 
replaced by multi-centered harmonic functions H°*(a#) = $; Sa + Gee 
satisfying the constraint 

H?) _ C2 _ G? 


= (89) 


u H2(2) CP GR 
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However, the process of determining kyng, and gyng in terms of charges and 
asymptotic moduli for a non-BPS multi-centered attractor is very different 
from its BPS counterpart. 

The reason is that the charge equations for a non-BPS multi-centered 
solution simplifies a great deal since 


[kn Bi, kNB, j] = 0 and [kn Bi, 9NB] =0 (90) 
guaranteed by the forms of (w”); and mAP. Therefore, the 5N equations 
decouple into N sets of 5 coupled equations: 


QnBi = So(knB,i) So (91) 


As in the single-centered non-BPS case, gyp does not enter the charge equa- 
tions (91), and its three degrees of freedom can be completely fixed by the 
given asymptotic moduli (xo, yo) and uo = 1 without using the zero NUT 
charge condition. More importantly, unlike BPS multi-centered solutions, the 
positions of centers x; do not appear in the charge equations (91), thus receive 
no constraint: all centers are free. Finally, since we are using the remaining 
3N +2 parameters {C°, z,u} to parameterize a N-centered attractor solu- 
tion under 5N constraints coming from charge equations (91), there need to 
be 2N — 2 constraints imposed on the D-brane charges. 

As in the BPS multi-centered attractors, solutions of C®® are given by the 
single-centered non-BPS solutions C°* = C°*(Q,z) with Q replaced by Qj. 
The solutions of twistors z and u are the same as the single-centered ones with 
charges Qng replaces by Q‘7',. Among the aforementioned 2N —2 constraints, 
N — 1 come from demanding that all centers have the same twistor z, which 
follows from the zero NUT charge condition at each center, and the other 
N — 1 come from demanding that they have the same twistor u as in (89). 
Solving these 2N — 2 constraints shows that all the charges {Qng} are the 


image of a single duality transformation Î ona multi-centered D4-D0 system 


{QNB 40,1}: M 
QNB, = T 'QNB,40,i (92) 
— The charges at different centers are all mutually local 
(QnB.i; QNB,j) = 0 (93) 


Like non-BPS single-centered attractors, the generic non-BPS multi-centered 
attractors cannot be generated via the “naive” harmonic function procedure, 
except for those with {QnB.i, h} being the image of a single Î on a pure 
D4-D0 system {QNB,40,i; hao}: 


ÎQNB aoi 
|z — zıl 


+ Îhao) (94) 


tna(x) = tġg(ÎQNB 10 > 5 


i 
In summary, the non-BPS multi-centered attractors are drastically dif- 
ferent from their BPS counterparts: there is no constraint imposed on the 
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positions of the centers, but instead on the allowed charges Qng: they have 
to be mutually local. The result is that the centers can move freely, and there 
is no intrinsic angular momentum in the system. 


6 Conclusion and Discussion 


In this talk, we summarized the construction of generic single-centered and 
multi-centered extremal black hole solutions in theories whose 3D moduli 
spaces are symmetric coset spaces. In this construction, all attractors, both 
BPS and non-BP%S, single-centered as well as multi-centered, are treated on an 
equal footing. The single-centered black hole attractors correspond to those 
null geodesics in M3p that are generated by exponentiating appropriate nilpo- 
tent elements in the coset algebra. The multi-centered black hole attractors are 
given by 3D solutions that live in certain null totally geodesic sub-manifolds 
of Msp. The construction of multi-centered attractors, even that of non-BPS 
ones, is merely a straightforward generalization of the single-centered con- 
struction. 

We presented a detailed computation in the theory of 4D M = 2 supergrav- 
ity coupled to one vector-multiplet, whose 3D moduli space is the symmetric 
coset space Gg(2)/SL(2,R)?. The attractor flow generators are third-degree 
nilpotent elements in the coset algebra. We explicitly constructed generic at- 
tractor solutions, both single-centered and multi-centered, and showed that 
while the BPS attractors can be generated from the attractor moduli via the 
“naive” harmonic function procedure, the generic non-BPS attractors cannot 
be generated this way. 

In the ny = 1 model, besides the BPS generator, there is only one extra 
third-degree nilpotent orbit to serve as non-BPS flow generators. Hence there 
is only one type of non-BPS single-centered attractor. In models with bigger 
symmetric moduli spaces, there should be more than one type of non-BPS 
generator. These would give rise to different types of non-BPS attractor flows, 
which might have different stability properties. 

All multi-centered non-BPS attractors constructed in this work follow from 
the ansatz in which 3D gravity is assumed to decouple from the moduli. The 
multi-centered non-BPS black holes are found to be very different from their 
BPS counterparts: the charges of all centers are constrained to be mutually 
local, while the positions of centers are completely free. Thus the non-BPS 
multi-centered attractor is not a “bound state” and carries no intrinsic angular 
momentum. 

We would like to construct true multi-centered non-BPS “bound states”, 
i.e., solutions with constraints on the positions of centers but not on the 
charges. There are two possible ways to achieve this. First, one could adopt 
a more general ansatz in which 3D gravity is coupled to the moduli. For 
axisymmetric configurations, the inverse scattering method could be used to 
perform an exact analysis. One could also search in models with bigger moduli 


Non-Supersymmetric Attractors in Symmetric Coset Spaces 25 


spaces. It is very likely that in bigger moduli spaces, there exist true multi- 
centered non-BPS “bound states” even within the ansatz with 3D gravity 
decoupled from moduli. We are also interested in the possibility of generating 
multi-centered non-BPS solutions with each center having different types of 
non-BPS generators ky p. 

Finally, with the hope of studying non-BPS extremal black holes in 4D 
N = 2 supergravity coupled to ny vector-multiplets with more generic pre- 
potential, we would like to generalize our method to non-symmetric homoge- 
neous spaces, and even to generic moduli spaces eventually. 
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